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We report a measurement of the difference (AAcp) between time-integrated CP-violating asym- 
metries in — >■ K" and — > tt+tt" decays reconstructed in the full data set of proton- 
antiproton collisions collected by the Collider Detector at Fermilab, corresponding to 9.7 fb~^ of 
integrated luminosity. The strong decay _D*+ — >• D'^tt^ is used to identify the charm meson at pro- 
duction as D° or . We measure AAcp = (-0.62±0.21 (stat) ±0.10 (syst))%, which differs from 
zero by 2.7 Gaussian standard deviations. This result supports similar evidence for CP violation in 
charm-quark decays obtained in proton-proton collisions. 
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The noninvariance of the laws of physics under the si- 
multaneous transformations of parity and charge conju- 
gation ( CP violation) is accommodated in the standard 
model (SM) through the presence of a single irreducible 
complex phase in the weak couplings of quarks. Generic 
extensions of the SM are expected to introduce additional 
sources of CP violation, which, if observed, could pro- 
vide indirect signs of new particles or interactions. To 
date, CP violation has been established in transitions 
of strange and bottom hadrons, with effects consistent 
with the SM interpretation Studies of CP viola- 

tion in charm decays provide a unique probe for non-SM 
physics. The neutral D system is complementary to the 
B and K sectors in that up-type quarks (electric charge 
+2/3) are involved in the initial state. Therefore, CP- 
violating effects probe the presence of down-type (charge 
— 1/3) new physics through charged-current couplings 
Q. However, Cf- violating effects are expected not to ex- 
ceed 0(10-2) in the SM fl, because charm transitions 
are well described by the physics of the first two quark 
generations. Indeed, no CP- violating effects have been 
firmly established yet in charm dynamics. 

Time-integrated CP-violating asymmetries of decays 
into CP eigenstates such as — !• tt+tt" and 
K^K^ probe non-SM physics contributions in the os- 
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dilation and penguin transition amplitudes. Oscillation 
indicates D'^-D transitions governed by the exchange 
of virtual heavy particles occurring before the D meson 
decay. Penguin decays are second-order transitions me- 
diated by an internal loop. Both amplitudes may be af- 
fected by the exchange of non-SM particles, which could 
enhance the size of the observed CP violation with re- 
spect to the SM expectation. In 2011, CDF reported 
CP- violating asymmetries compatible with zero within a 
few 10""^ uncertainty in these decays, along with a mea- 
surement of the difference AAcp of CP- violating asym- 
metries in K^K^ and tt+tt", also con- 
sistent with zero |8|]. Shortly after, LHCb reported a 
more precise determination of the difference, which is 
3.5(7 different from zero Q. If established, this result pro- 
vides the first evidence for CP violation in charm dynam- 
ics, with a size larger than most SM expectations [lo| . 
Among the quantities sensitive to CP violation in charm 
mesons, /^Acp can be measured with good accuracy be- 
cause many common systematic uncertainties cancel. In 
addition, lS.Acp could be maximally sensitive to CP vi- 
olation since the individual asymmetries are expected to 
have opposite signs, if the invariance of the dynamics 
under interchange of d with s quarks is approximately 
vahd 0. 

In this Letter, we report a measurement of the differ- 
ence of CP-violating asymmetries in D° — >■ K'^K^ and 
— !■ TT+TT" decays reconstructed in the full data set 
of 1.96 TeV proton- antiproton collisions collected by the 
Collider Detector at Fermilab, corresponding to 9.7 fb~^ 
of integrated luminosity. In addition to the increase in 
the size of the data set with respect to the previous mea- 
surement , we optimize the selection for a measurement 
of the difference of asymmetries, reaching a sensitivity 
competitive with the LHCb result . 

For each decay mode, the CP-violating time- 
integrated asymmetry between decays of states produced 
as 1?° and I?" is defined as 



N{D^ h+h-;t) - N{D^ h+h-:t) 
N{DO h+h-;t) + N{W h+h-:ty 



(1) 



where h identifies a charged kaon or pion. The asym- 
metry can receive contributions from any difference in 
partial decay widths between and mesons (direct 
CP violation) and both the difference in mixing proba- 
bilities between and D mesons and the interference 
between mixed and unmixed decays (indirect CP viola- 
tion). Due to the slow mixing rate of charm mesons, the 
asymmetry is approximated to first order as the sum of 
two terms [sl, 



{t) 



Acp{h+h~) = A^^(/i+/i-) + ^ A'^f,{h+h-). 

T 



(2) 



The first term arises from direct CP violation and de- 
pends on the decay mode, the second from indirect CP vi- 
olation and is nearly independent of the decay mode 11 



The average decay time of the sample used in the mea- 
surement {t) depends on the detector acceptance as a 
function of decay time and r is the D° lifetime. To 
measure each individual asymmetry, we determine the 
number of detected decays of D° and D and use the 
fact that primary charm and anticharm mesons are pro- 
duced in equal numbers by the CP-conserving strong 
pp interactions. We require the D candidate to be pro- 
duced in the decay of an identified D*'^ or D*~ meson 
to determine whether the decaying state was initially 
produced as a or a D° meson. Flavor conserva- 
tion in the strong-interaction decay of the D*^ meson, 
D*+ D°{^ h+h-)7r+ and D*- d"{^ h+h-)Tr^, 
allows identification of the initial charm fiavor through 
the sign of the charge of the low-momentum tt meson 
(soft pion, TTs). The observed asymmetry, A{h^h^) — 
Nohs{D° -)■ h+h-) - Nohs(D° h+h-) /[anm], is the 
combination of the contributions from CP violation and 
from the detection asymmetry between _D*+ and D*^ 
mesons, due to different reconstruction efficiency for posi- 
tive and negative soft pions. The combination is linear for 
small asymmetries, yl(ft,+ /i^) = Ac'p(/i^/i^)+i5(7rs)'' . 
The instrumental asymmetry is due to differences in in- 
teraction cross sections with matter between positive and 
negative low-momentum pions and the geometry of the 
CDF tracking system (see below). The combined effect 
of a few percent cancels in the difference of asymmetries 
between K^K^ and tt+tt^ decays, 

AAcp =A{K+K-) - yl(7r+7r") 

=Acp{K+ K-) - Acp{^+^-) 



--AAtp^ 



■^'cpA(i)/r. 



Kinematic differences between the K^K~ and tt+tt" 
decays result in a fractional 10% difference in average 
decay-time of the sample, A{t)/T = 0.27 ± 0.01, mea- 
sured through a fit to the decay-time distribution of 
background-subtracted signal candidates Therefore, 
most of the indirect CP-violating asymmetry cancels in 
the subtraction and AAcp approximates the difference 
in direct CP- violating asymmetries of the two decays. 

The CDF II detector is a multipurpose magnetic spec- 
trometer surrounded by calorimeters and muon detec- 
tors. The detector components relevant for this analysis 
are briefly outlined below; a more detailed description is 
in Ref. '&\. A silicon microstrip vertex detector and a 
cylindrical drift chamber immersed in a 1.4 T axial mag- 
netic field allow reconstruction of charged-particle tra- 
jectories (tracks) in the pseudorapidity range |?7| < 1.0. 
The vertex detector contains seven concentric layers of 
single- and double-sided silicon sensors at radii between 
1.5 and 22 cm, each providing a measurement with up 
to 15 (70) /im resolution in the cj) (z) direction [121 . The 
drift chamber has 96 measurement layers, between 40 
and 137 cm in radius, organized into alternating axial 
and ±2° stereo superlayers [l3|. A 35° tilt angle between 
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the drift chamber cell orientation and the radial direction 
facilitates track finding, but induces charge-dependent 
detection asymmetries of up to a few percent for low- 
momentum charged particles [1, ■ The component of 
a charged particle's momentum transverse to the beam 
(Pt) is determined with a resolution of crpj./p^ « 0.07% 
{pt in GeV/c), corresponding to a typical mass resolu- 
tion of 8 MeV/c^ for a two-body charm- meson decay. 

The data are collected by a three-level online selection 
(trigger) system. At level 1, tracks are reconstructed 
in the transverse plane. Two oppositely-charged parti- 
cles are required, with reconstructed transverse momenta 
Pti, PT2 > 2 GeV/c, the scalar sum pTi + Pt2 > 5.5 
GeV/c typically, and an azimuthal opening angle Acf) < 



90° [15|. At level 2, tracks are combined with silicon 



hits and their impact parameter d (transverse distance 
of closest approach to the beam line) is determined with 
45 /zm resolution (including the beam spread) and typi- 
cally required to be 0.12 < d < 1.0 mm A slightly 
tighter opening-angle requirement of 2° < A0 < 90° is 
also applied. Each track pair is then used to form a neu- 
tral D candidate, whose flight distance in the transverse 
plane projected onto the transverse momentum L^y is 
required to exceed 200 /im. At level 3, the selection is 
confirmed on events that are fully reconstructed by an 
array of processors. 

The offline selection is based on a more accurate de- 
termination of the same quantities used in the trigger 
and is detailed in Ref. [1, [l3]. Here we only describe 
the improvements specifically aimed at enhancing the 
sensitivity to AAcp- By exploiting the highly accurate 
cancellation of instrumental effects in AAcp, the selec- 
tion of Ref. is loosened to significantly increase sig- 
nal efficiency. Requirements on the minimum number of 
measurement points for reconstructing tracks are loos- 
ened, and the transverse momentum threshold for D de- 
cay products is lowered from 2.2 to 2.0 GeV/c, nearly 
doubling the yield. Asymmetries from charm mesons 
produced in B meson decays (secondary charm) may in- 
troduce a common bias in the individual CP-violating 
asymmetries, which cancels in the difference of asymme- 
tries. Secondary charm decays are therefore not excluded 
from the analysis, providing a 12% additional signal. An 
additional 10% of events are contributed by secondary 
trigger selections that were not used in Ref. Finally, 
the additional integrated luminosity contributes approx- 
imately 20% more events. As a result, the expected av- 
erage resolution on the difference of asymmetries is im- 
proved by approximately 50% with respect to Ref. [§]. 

The reconstruction of signal candidates is entirely 
based on tracking information, with no use of particle 
identification. Two tracks from oppositely-charged par- 
ticles, with pion or kaon assignment, are fit to a com- 
mon decay vertex to form a neutral D candidate. A 
low-momentum (px > 400 MeV/c) charged particle is 
associated with each D candidate to form the charged 



D* candidates. This allows identification of the charm 
meson at production as a D'^ or a D and strongly re- 
jects background, albeit with an 85% reduction in sig- 
nal yield. In the few percent of cases in which multiple 
candidates per event are reconstructed, one is randomly 
chosen for further analysis. The resulting sample con- 
tains approximately 591 000 — !■ K^K^ candidates, 
619 000 if K+K- candidates, 270 000 D° tt+tt" 
candidates, and 279 000 D tt+tt" candidates. Many 
kinematic distributions are compared for the K^K^ and 
TT+TT^ samples. Small differences are observed in the 
distributions of the soft pion's impact parameter, trans- 
verse momentum, and pseudorapidity. The final candi- 
dates are reweighted to make these distributions equal, 
by using as a weight a smooth function extracted from 
a fit to the ratio of K^K^ and tt+tt^ distributions. 
The data show that the weight function factorizes into 
the product of three one-dimensional functions and typi- 
cally does not exceed 10%. The reweighting ensures that 
any kinematics-dependent instrumental asymmetry can- 
cels in the difference of observed asymmetries. 

The observed asymmetry in each sample is determined 
from a simultaneous fit of the D^ir^ and D Tr~ binned 
mass distributions of candidates restricted to the signal D 
region, defined as those with mass within 24 MeV/c^ {3a) 
of the known mass. The D'^tt mass is calculated using 
the vector sum of the momenta of the three particles to 
determine the D* momentum and the known D'^ mass. 
This quantity has the same resolution advantages of the 
more customary M{h~^h~Tr) — M{h'^h~) mass difference, 
but it is independent of the mass assigned to the decay 
products. 

The 7T~^Tr~ sample is dominated by the signal of D*- 
tagged decays, a background of real decays as- 
sociated with random pious or random combinations of 
three tracks (combinatorics), and a 0.93% contamina- 
tion of the high-mass tail of the K~7t~^ signal 
misreconstructed as a tt+tt^ final state. In the K^K^ 
sample, an additional background is contributed by mis- 
reconstructed multi-body charm meson decays, domi- 
nated by the D° /i^tt+tt" and the 13° h^l+vg con- 
tributions, where ^ is a muon or an electron. The func- 
tional form of the signal mass shapes is determined from 
simulation, with parameters tuned in a low-background 
sample of 12.5 x 10^ L»*"tagged D° K-TT+ decays [1]. 
The data indicate a small asymmetry between signal 
shapes of Z3*+ and D*^ decays, which can be attributed 
to differences in the tracking resolution for positive and 
negative soft pions. This asymmetry is included in our 
fit model [1]. The shape of the combinatorial compo- 
nent is obtained and fixed from data, by forming arti- 
ficial D* candidates where each candidate is associ- 
ated with soft pions of all candidates found in different, 
randomly chosen events. The functional form of the mis- 
reconstructed decays is extracted from samples of simu- 
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FIG. 1: Distributions of _D°7r^ mass with fit results overiaid for (a) D" — >■ n'^Tv decays, (b) D° — !> n^n decays, (c) 
D° K+K- decays, and (d) if K+ K~ decays. 



lated inclusive charm meson decays, and its parameters 
are fit together with the desired asymmetries 0. The 
Ktt tail is not included in the tttt fit, but is accounted 
for in the systematic uncertainties. In the tt+tt" sample, 
the parameters determined by the fit are the asymme- 
try between _D*+ and D*^ yields and the relative sizes 
between the signals and the combinatoric background 
components. In the K~^K~ sample, the fit also deter- 
mines the relative sizes and values of the shape parame- 
ters of the mis-reconstructed component. The fit allows 
for asymmetries between the numbers of combinatorial 
and mis-reconstructed background events in the D*^ and 
D*~ samples. For each final state, we minimize the total 
for the D*~^ and D*~ samples and obtain the results 
shown in Fig. [T] The fits show agreement with data, and 
the observed asymmetries are 

A{t:+tt-) = (-1.71 ± 0.15 (stat))%, (3) 
A{K+K-) = (-2.33 ± 0.14 (stat))%. (4) 

Both asymmetries are dominated by the detector-induced 
contribution. Fits including extreme variations of the 



signal and background models yield significantly larger 
values of reduced with minimal variations in the ob- 
served asymmetries. As a consistency check, we compare 
the results of the measurement obtained in independent 
subsamples chosen according to the soft pion's direction 
in the four quadrants of the tracking volume, different 
data-taking periods (early and late data), or splitting 
the present sample in the subsample of events used in 
Ref. [§] and the complementary sample, in which we ob- 
serve AAcp = -0.74 ± 0.27. The results show a high 
level of consistency, with reduced between observed 
asymmetry differences of 4.4/3, 0.38/1, and 0.46/1, re- 
spectively. 

Most systematic effects cancel in the subtraction of 
asymmetries. Residual higher-order instrumental effects 
that do not cancel are estimated to contribute less than 
±0.009% to AAcp, based on simulations in which known 
instrumental asymmetries are varied as functions of the 
kinematic variables. The impact of possible residual mis- 
modeling of the mass shapes used in fits is evaluated 
by repeating the measurement using extreme variations 



7 



of the model, as derived from data, and contributes to 
AAcp by less than ±0.020%. A dominant systematic 
uncertainty of 0.1% arises from the possibility that signal 
and background shapes differ between the D*+ and D*~ 
samples. This effect is assessed by repeating the fit on 
data using various modifications of the fit shapes in which 
independent parameters are used for _D*+ and D*~ sam- 
ples. The effect of the Ktt tail in the tttt signal induces 
a systematic uncertainty of 0.013% that is the product 
of the size of the contamination (0.93%) times the 3% 
observed asymmetry of the — J> K~Tr'^ decay. The im- 
pact of the statistical uncertainties associated with the 
kinematic reweighting is negligible. 
The final result, 

AAcp = (-0.62 ± 0.21 (stat) ± 0.10 (syst))%, (5) 

is consistent with and supersedes the previous CDF 
determination of AAcp = (-0.46 ± 0.31 (stat) ± 
0.12 (syst))% By adding in quadrature the un- 

certainties, assumed to be independent and Gaussian- 
distributed, the difference of asymmetries deviates from 
zero by 2.7 standard deviations, strongly indicating the 
presence of CP violation in the decays of D'^ mesons. 
This result is consistent with the LHCb measurement 
obtained in pp collisions, AAcp — (—0.82 ±0.21 (stat) ± 
0.11 (syst))% i, with comparable accuracy and less 
than la difference in central value. The combined results 
of the two experiments provide substantial evidence for 
CP violation in the charm sector with a size larger than 
most predictions [l^ , possibly suggesting the presence of 
non-SM dynamics. More precise determinations of the in- 
dividual asymmetries in Z)" — ;> tt+tt" and — J> K^K~ 
decays and extension of the precise experimental explo- 
ration to other charm decays may help in understanding 
whether the observed effect can be attributed to signifi- 
cant hadronic corrections to the SM weak amplitudes or 
to new, non-SM sources of CP violation [l3| . 
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